Optical diffraction spectroscopy of excitons in uniaxially-strained GaN films 
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The degenerate four- wave mixing spectroscopy of uniaxially strained GaN layers is demonstrated 
using colinearly polarized laser pulses. The nonlinear response of FWM signal on exciton oscilla- 
tor strength enhances the sensitivity for polarized exciton, allowing for mapping out the in-plane 
anisotropy of the strain field. The observed high-contrast spectral polarization clearly shows fine 
structure splittings of excitons, which are also confirmed in the change of quantum beating periods 
of time. 



Influence of strain-fields on the optical properties of 
semiconductors has been studied extensively for many 
years. In particular, the presence of uniaxial strain gen- 
erates the polarization and splitting of excitons, which 
provide useful information on exciton fine structures 
(EFSs)±i2iM. However, it is only recently that the EFS 
in wurtzite GaN was observed experimentally, where the 
exchange splitting was measured by M. Julier et at. as 
0.6±0.1 meV using reflectivity measurements on a uniax- 
ially strained GaN epilayer 5 , and by P. Paskov et al. as 
0.58±0.05 meV using photoluminescence (PL) measure- 
ments on thick GaN 6,7 . The exchange interaction con- 
stant of 7 evaluated by a Hamiltonian based on quasi- 
cubic approximation shows a good correspondence in 
both measurements and is consistent with the exponen- 
tial dependence on interatomic distance for III-V com- 
pounds. However, since the EFS is usually comparable 
to the exciton linewidth, more precise measurements are 
needed for detailed analysis. The use of nonlinear spec- 
troscopy often provides significant enhancement of the 
signal sensitivity, which is advantageous for investigating 
a small amount of spectral changes such as those of EFSs. 

In this study, we report on the polarization-resolved 
mapping of exciton lines in a uniaxially strained GaN 
layer using degenerate four-wave mixing (FWM) spec- 
troscopy. Since the diffraction signal in FWM increases 
as the fourth power of the exciton oscillator strength, the 
exciton spectra show strong polarizations resulting from 
the uniaxial strain in the sample. The highly polarized 
spectra allow us to determine the precise splitting en- 
ergy between excitons, which shows a good agreement 
with that obtained from the oscillation period of quan- 
tum beats. 

The two-pulse FWM experiments in reflection geom- 
etry were performed using a frequency doubled, mode- 
locked tunable Tiisapphire laser with the spectral width 
of 13 meV (FWHM) 8 . The excitation pulses with the 
same intensities were colinearly polarized and superim- 
posed onto a sample surface using a lens (f= 200 mm). 
The FWM signal in the 2k 2 -k\ direction was spatially 
selected by an iris, and was spectrally resolved by a 30- 
cm-long monochromator with a CCD detector. The de- 



lay time Ti2 between the pulses is defined as positive 
when the k\ pulse precedes the k 2 pulse. The sample is 
mounted in a closed-cycle helium cryostat and all mea- 
surements were carried out at 10 K. 

To understand the enhancement of sensitivity for ex- 
citon polarizations, we consider the FWM signal in 
a two-level system using the phase-space filling (PSF) 
nonlinearit}iifl. For simplicity, we neglect additional 
many-body Coulomb interactions. Within the framework 
of the optical Bloch equation for the ^-function pulses 
with the linear polarizations E^ 2 and Ek ± , the third-order 
optical polarization _ k at rn is described b}A 

p il- kl « -\®{t-T 12 )Q{T 12 )ElE*M^ q , 

where \i q is the dipole matrix elements of linearly polar- 
ized excitons. As a consequence, the FWM signal Ifwm 
with respect to the oscillator strength \ji q \ 2 is expressed 
by 

Ifwm(t 12 ) (x J dt\P^l_ ki -/i*| 2 , 

indicating the fourth power dependence of |/i g | 2 . The 
FWM signals thus provide a sensitive probe for linearly 
polarized excitons. 

The sample used in this study is a high-quality 2.3 fim- 
thick GaN layer on an (11-20) a-plane sapphire substrate 
grown by metal-organic chemical vapor deposition^!. 
The details of the sample growth conditions were pre- 
sented in Ref. 11. In the sample, built-in uniaxial strain 
exists due to the anisotropic thermal expansion coeffi- 
cients of the a-plane sapphire. Figure 1(a) shows the ex- 
perimental geometry with respect to the crystallographic 
axis of the sample. The polarization angle of the inci- 
dent pulses, which can be changed by rotating the half- 
wave plate (HWP) in front of the focusing lens, is mea- 
sured from the [11-20] axis parallel to the direction of the 
uniaxial stress. Figure 1(b) illustrates an overview of the 
optically allowed states of A-exciton (X A ) and B-exciton 
(X B ) when the laser light is polarized perpendicular to 
the growth axis. In normal wurtzite GaN, there exists 
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twofold degenerate spin-singlet states. The reduction of 
symmetry due to the uniaxial strain lifts the degeneracy 
and polarizes the optical transitions allowed in x and y, 
where x and y are along the [11-20] and [1-100] axes, 
respectively. 
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FIG. 1: (a) Experimental geometry. The colinear polariza- 
tion angle of the incident pulses is defined from the [1-100] 
axis of GaN. (b) Schematic illustration of F$ exciton levels in 
wurtzite GaN. Uniaxial in-plane strain reduces the symme- 
try from Cq v (hexagonal) to Civ (orthorhombic). (c) Typical 
FWM spectra of GaN layer on a-plane sapphire (top), GaN 
layer on c-plane sapphire (middle), and free-standing GaN 
(bottom), at r±2 = s. In the top spectra, the polarization 
angle is set to be tt/4, and the dashed line indicates the 
excitation spectrum. 



analysis of A. Shikanai et. a! 13 , we obtained A cr = 30 
meV and A so = 13 meV. 
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FIG. 2: (color). Contour plot of the FWM spectra for a GaN 
on a-plane sapphire taken by rotating the polarizaion at m 
— 0. The spectra consists of peaks centered at approximately 
3508.6 meV (X A ) and 3499.0 meV (X s ). The polarization 
angle of incident pulses was rotated in 1-degree steps and the 
integration time for collecting the signal at each was 1 sec. 



TABLE I: Exciton energies obtained from FWM spectra at 
10 K. 



Sample 


X A (eV) 


X B (eV) 


X B -X A (meV) 


Bulk GaN 


3.4816 


3.4871 


5.5 


GaN on c-sapphire 


3.4912 


3.4987 


7.5 


GaN on a-sapphire 


3.4990 


3.5086 


9.6 



We first consider the biaxial strain in order to deter- 
mine several important parameters of crystal-field split- 
ting A cr and spin-orbit splitting A so . Figure 1 (c) shows 
typical FWM spectra at t\2 = 0. For comparison, FWM 
spectra of a 70 /im-thick free-standing GaN and a 2.3 
/xm-thick GaN layer on a (0001) c-plane sapphire are 
displayed in the figure. In the spectra, the lower- and 
upper-energy peaks are X A and X s , respectively, each of 
which shows the systematic energy shifts depending on 
the sample. The exciton energies and their separations 
obtained from the Lorentzian fitting to the spectrum in 
each sample are listed in Table I. As has been reported 
by several group a 13 i 14 , the observed shifts reflect the bi- 
axial strain associated with the mismatch of the lattice 
constants and of the thermal expansion coefficients be- 
tween GaN and substrates in each sample. Following the 



The ability to map out the uniaxial strain-field is high- 
lighted in Fig. 2, where we show a contour-plot of FWM 
spectra consisting of X A and X s , recorded at the po- 
larization angle of the incident pulses with a 1-degree 
step. The figure shows a striking feature of the anti- 
corr elated polarizations in each exciton line, signifying 
the reduced symmetry due to the uniaxial strain. This 
optical anisotropy has been reported previously using 
the reflectance and the PL measurements^*!. Note that 
compared with the case of linear spectroscopy, there is 
a considerable magnitude of polarization in which the 
nonlinear response in the FWM measurement increases 
by a factor of ~10, allowing us to determine the uniax- 
ial strain-field without fitting. Besides the nonlinearly 
enhanced exciton polarization, the contour-plot clearly 
shows the energy shifts originating from the EFS of two 
linearly polarized excitons. The lower-energy component 
of X A and the upper-energy component of X B are active 
at = 0, 7r corresponding to the [1-100] axis (y-direction 
in Fig. 1 (b)), whereas the upper-energy component of 
X A and the lower-energy component of X B are active at 6 
= 7r/2, 37r/2 corresponding to the [11-20] axis (x in Fig. 
1 (b)). Thus, the polarization spectra provides a type 
of crystalline analysis achieved in X-ray diffraction spec- 
troscopy. Contrary to the X-ray diffraction, the FWM 
analyzes the sample using safe and simple experimental 
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FIG. 3: Plots of (a) normalized peak intensities and (b) en- 
ergies on A-exciton (red) and B-exciton (blue) lines obtained 
from the Lorentzian fitting to the FWM spectra at tyi =0.2 
ps for various polarizations. Calculated polarization depen- 
dences of (c) FWM intensities and (d) peak energies for each 
exciton line. 



apparatus. We stress that this technique can be applied 
extensively to anisotropically polarized exciton systems 
with uniaxial strain. 



TABLE II: Fitting parameters for polarized FWM. E a and 
(a=A,B) are exciton transition energy and linewidth 
(FWHM), respectively. 



Polarization 


E A (eV) 


T A (meV) 


E B (eV) 


T B (meV) 


X 


3.4993 


2.1 


3.5082 


1.8 


Y 


3.4987 


2.3 


3.5091 


2.5 



For qualitative analysis with great precision, we made 
to fit the spectra with a combination of two Lorentzian 
functions. A Lorentzian lineshape is predicted in the case 
of a homogeneously broadened system, which will be con- 
firmed in the time-evolutions of the signal. Figure 3 (a) 
and (b) plot the peak intensities and the energies for each 
exciton line at various 0. Note that we evaluate the data 
obtained at t\2 ~ 0.2 ps in order to reduce the contri- 
butions from efficient many-body effects at t\i ~ 0. In 
Fig. 3 (a), the variations of the FWM intensities deviate 
from the simple sinusoidal functions but agree well with 
the fourth power of them [sin 4 (#)], providing proof that 
the PSF nonlinearity accounts for our detection. The 
best fit to the data presented in Fig. 3 (c) gives the 
ratios of polarized FWM intensities for X A and X B as 



jAxjjAy = 19.0^(2. 1) 4 and I B v / I Bx = 17.2^(2.0) 4 , re- 
spectively. This results in a trapezoidal variation of the 
spectral peak energies (Fig. 3(b)). Based on the ratios of 
polarized FWM intensities, we can reproduce such polar- 
ization dependences as shown in Fig. 3 (d). The param- 
eters obtained from these analysis are listed in Table II. 
The splitting energies between the colinearly polarized 
EFS in each exciton are E B -E A = 8.9 meV (x) and 10.4 
meV (y). The average energy shift of the EFS is thus 
~0.8 meV. Using the specific parameters obtained from 
the biaxial strain, we are now ready to evaluate the ex- 
change interaction constant 7 for this sample. Following 
the analysis presented by M. Julier et a/. 5 , we obtain 7 
= 0.6 meV using the uniaxial strain energy 63 = —1.5 
meV. 
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FIG. 4: (a) Contour plot of FWM spectra at = tt/4 for 
various T12. FWM intensities for (b) X# and (c) Xa as func- 
tions of delay T12 at different incident polarizations 0. The 
gray lines are guides for the eye, indicating the changes in 
quantum beating oscillation Q for various 0. 



In addition to the spectral evolutions, FWM exper- 
iments allow us to examine the temporal dynamics of 
excitons in EFS. Figure 4 (a) shows a contour plot of 
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FWM as a function of T12, at 6 =7r/4. As is well known, 
the simultaneously excited two excitons result in a quan- 
tum beating oscillation, of which the oscillation period Vt 
corresponds to the energy separation between excitons. 
In Figs. 4 (b) and (c), we plot the peak intensities of 
each resonance at = 0,7r/4,7r/2. The polarization de- 
pendence of Q AjB (0) in each polarization configuration is 
clearly seen, and is, of course, comparable in each exci- 
ton line (Q A (6) « Q B (0)). Assuming the homogeneously 
broadened excitons, the FWM signal at T12 can be ex- 
pressed as 

[1 + a A ' B (0) cos(Q(0)r 12 - ^ B )} exp(-2r 12 /T 2 A ' B (#)). 
In this fitting function, a A,B denotes the beating contrast 

A B 

and T 2 ' denotes the dephasing time, both of which de- 
pend on 6. Neglecting the many-body Coulomb inter- 
actions, we consider the quantum beats dominated by 
the two pairs of colinearly polarized excitons, X. Ax -~X. Bx 
and X Ay -X By . In this assumption, the results obtained 
at = and tt/2 reflect the natures of X. x and re- 
spectively. The exact phase (j) AjB is slightly different in 
each polarization mainly due to the Coulomb interac- 
tions between excitons-^. For simplicity, we also neglect 
these contributions. From the fitting to the data, we ob- 
tain Q x = 13.6 rad./ps and Q y = 15.7 rad./ps. These 
parameters correspond to the splitting energies between 
the colinearly-polarized EFS in each exciton of E B -E A = 
8.95 meV (x) and 10.3 meV (y), coinciding with those 
obtained from the polarization-resolved spectra. 

The dephasing times in each exciton are T Ax = 0.58 ps, 
T$ v = 0.60 ps, T Bx = 0.70 ps, and T^ y = 0.58 ps, which 
are comparable with those obtained from the Lorentzian 
fits to the FWM spectra. Thus, we conclude that exciton 
lines are homogeneously broadened in our sample. 

Also, the beating contrasts a A,B (0) are a A (0) w 
a B (0) « 0.6 and a A {ii/2) « a B {it/2) « 0.9. The po- 
larization dependence is mainly due to the difference in 
relative oscillator strength between the excitons. On the 
basis of a noninteracting three-level system, the beating 
contrast is given by 16 

1 ; |/m(0)I 4 + IM0)I 4 ' 



According to the spectral evolutions, we assume the 
relative exciton oscillator strengths to be |^^i(0) | 2 = 
0.31, \fi B (0)\ 2 = 1.0, |/MW2)| 2 = 0.65, and |^(tt/2)| 2 
=0.49. As a result, we obtain the beating contrasts as 
a A ' B (0) « 0.57 and a A ^ B (ir/2) « 0.96, which give ex- 
cellent agreement with those evaluated using the time 
evolutions. 

Finally, We comment on the splitting of EFS in each 
exciton. In Fig. 3(b), the shifts of EFS energy are 

E Ax _ E Ay = q g me y and E By _ E Bx = q g meV? ghow _ 

ing a difference in splitting energy between X A and X B . 
This difference between excitons is consistently observed 
throughout all our FWM experiments, and becomes even 
more apparent in Fig. 2, where, in the absence of any fit- 
ting, one can see a clear difference in the distributions of 
each exciton line. Furthermore, we have obtained a sim- 
ilar difference in the reflectivity measurements, exclud- 
ing the contributions of many-body effects. Therefore, 
we conclude that the different splittings originate from 
the difference in EFS between excitons of our sample, 
suggesting the breakdown of the quasi-qubic approxima- 
tion in wurtzite GaN 17 . Further investigations including 
theoretical supports are needed to clarify the difference. 
However, we emphasize the importance of the enhanced 
sensitivity for exciton polarizations achieved in this mea- 
surement that allows us to account for the small amount 
of spectral changes in EFSs. 

In conclusion, a combination of the fourth power de- 
pendence of exciton oscillator strength on a four-wave 
mixing signal and the built-in uniaxial strain of the sam- 
ple allow not only for a type of crystalline analysis by X- 
ray diffraction spectroscopy, but also for studying exciton 
fine structures and their temporal dynamics. This tech- 
nique can be applied to optical characterization for po- 
larized exciton systems with an appropriate anisotropic 
strain field. 
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